Stability of hydrogenated group-IV nanostructures: magic structures of diamond 

nano crystals and Silicon quantum dots 
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We have developed an effective model to investigate the energetic stability of hydrogenated group- 
IV nanostructures, followed by validations from first-principles calculations. It is found that the 
Hamiltonian of X m H n (X=C, Si, Ge and Sn) can be expressed analytically by a linear combination of 
the atom numbers (m, n), indicating a dominating contribution of X— X and X— H local interactions. 
As a result, we explain the stable nanostructures observed experimentally, and provide a reliable 
and efficient technique of searching the magic structures of diamond nanocrystals(Dia-NCs) and 
Silicon quantum dots(SiQDs). 

PACS numbers: 61.46.-w, 61.46.Df, 68.65.-k 



Semiconductor nanocrystals have been greatly attrac- 
tive and intensively investigated p], Q. These nanomate- 
rials extend the physics of reduced dimensions and offer 
the opportunity for fundamental study of the regime be- 
tween nanostructure and bulk states [3j, which have also 
brought wide applications as nanoscale electronic and op- 
tical devices fluorescent biological labels [5], quantum 
computation media [6] etc. Besides II- VI and III-V semi- 
conductor compounds 0, various hydrogenated group-IV 
(C[|,[9|, Si^El, and Ge[H) nanocrystals with sp 3 hy- 
bridizations have been synthesized and isolated, whose 
optical response depends on particle size, shape and sym- 
metry. Theoretical studies |13l-[l5|, have focused on opti- 
cal properties of these nanocrystals and employed vari- 
ous methods for accurate calculations of adsorption spec- 
trum. 

The determination of stable structures, which dom- 
inate the optical properties of nanocrystals [H @, is not 
well understood so far. Previous calculations constructed 
Dia-NCs[l5| according to the synthesis experiments, 
which indicated that larger members of the series have 
smaller surface-to- volume ratios and lower hydrogen-to- 
carbon ratios [8|. According to the Wulff energy, the poly- 
crystalline wire of five-fold symmetry is more stable than 
single-crystal types [16], for silicon nanowires with the di- 
ameter less than 6 nm. Recent studies [17], Gil searched 
the magic structures of silicon nanowires using genetic 
algorithm, in which the energies are calculated with clas- 
sical potential in Hansel- Vogel (HV) formulism to reduce 
the time cost. 

It is challenging to determine the stable configuration 
from numerous possible candidates and there are three 
main obstacles for searching magic structures of X m H n 
(X=C, Si, Ge and Sn): 1) the accurate calculation of 
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the total energy is necessary but often computationally 
expensive; 2) over many isomeric structures should be 
considered and the number of these structures increases 
sharply as m and n increase; 3) the minimum of n is 
not clear for a certain m, though the maximum of n is 
trivially 2m + 2. 

In this letter, we investigate hydrogenated group-IV 
nanocrystals by both model analysis and first-principles 
approaches. We proposed an effective model and gave 
an analytical expression of the Hamiltonian for X m H n 
(X=C, Si, Ge and Sn) with the numbers of atoms (m, n), 
as is confirmed by the first-principles calculations. Our 
finding provides an efficient and reliable avenue of search- 
ing magic structures of Dia-NCs and SiQDs, which ex- 
tends our understanding on experiment observed stable 
nanocrystals. 

In our model, we assume the Hamiltonian of X m H n 
(X=C, Si, Ge and Sn) to be as 



H = 2^ (Hint + Ho) - rn/ix - n/i H (1) 



where Hint an d Ho are the contributions from interac- 
tions and self energies, and fixi/^B.) is the chemical po- 
tential for X(H) atom. The sum runs over all the group- 
IV atoms and we have Ho = Mx +_PiMH ? where pi is 
the number of H atoms in saturated group of the ith 
X atom and ^x (^u ) is the isolated atomic energy for 
X(H) atom. Hint includes the energy contributions from 
bonded X-atom pair(— Ex-x) and the saturated group 
of the ith X atom(— E-xn p .)- As shown in the inset of 
Fig.l, every X atom has four nearest neighbors and ev- 
ery X— X bond is shared by two X atoms. This leads the 
energy contribution corresponding to the ith X atom to 
be —2 Ex-x f° r Pi — 0, an d —1.5 Ex-x — ^-xh for Pi = 1 
analogically. Thus, H[ nt = -(2 - 0.5pi)E X -x ~ E-xu Pi • 

We assume that the interaction between X and H 
atoms is localized and thus E xr v . = PiEx-n, Pi = 




FIG. 1: (Color online) Possible (m, n) for X m H n . The convex 
corresponds to the local minimum and magic structures of 
hydrogenated nanostructures. 



0,1,2,3. With n = E™iP;> we found that Eq. (1) 
can be written as H = am + bn, where a = (—2Ex~x — 
Mx + MXo) and 6= (-E X -k + 0.5£ X -x ~m + Mh )- Tnis 
implies that the energy of X m H n is determined by the 
interaction strength of X— X and X— H, chemical poten- 
tials, and the atom numbers (m,n). 

Instead of scanning; the parameter space (a, b)[19], the 
Simplex method pol-[22j is efficient to determine the com- 
binations of lowest energy from a set of possible inte- 
ger combinations (m,n), which indicates that the stable 
configurations correspond to corners, edges and faces [2 2\. 
We will obtain the ground states if we construct enough 
restricting inequalities and solve the corresponding lin- 
ear euquations[22]. In our case, we have m > 1 and 
n < 2m + 2. The key task is to determine the lower 
limit of n for a certain m, for which it is difficult to find 
out the expected inequalities. We search the least n as 
follows (23j: i ) starting from one of the stochastic config- 
uration of XmHn (n < 2m + 2), exchange X and H atoms 
and saturate the configuration with H atoms when nec- 
essary; ii ) we accept the new configuration when the 
H atom number is non-increasing, otherwise the new ge- 
ometry will be accepted with the probability of l/\dn\, 
where dn is the increment of H-atom number. We find 
that n will converge into the minimum after hundreds of 
iterations. 

Figure 1 shows the upper and lower limit of n as a 
function of m. It is not a standard Simplex because both 
m and n will be increasing with an increasing nanocrystal 
size of X m H n . However, we will obtain a convex quad- 
rangle ABCD (colored in red) if a restriction of m < 8 
is considered. The convex A (1, 4) and B (6, 12) are 
corresponding to stable configurations, while C (8, 16) 
and D (8, 18) are not since they are induced by the ar- 
tificial restriction of m < 8. Thus, XH4 and X6H12 will 
be stable configurations for group-IV nanocrystals. Ana- 
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FIG. 2: (Color online) Cohesive energies of Dia-NCs and 
SiQDs as a function of the H/(C, Si) ratio. Inset shows the 
linear dependence between the reaction heat of Dia-NCs and 
hydrogen numbers. 



logically, we will obtain a new convex quadrangle (color 
in red and yellow) and find another stable configuration 
of X10H16 in place of X6H12, if we consider a restriction 
of m < 12. Besides, we find that X14H20, X18H24 and 
X22H28 are also stable configurations. It should be noted 
that, except for XH4, all other nanocrystals are meta- 
stable states, since they are local convex ascribed to the 
size confinement. As is known, the size of nanocrystals 
increases with increasing reaction time as more material 
is added to the surfaces. 

Till now, we have settled the main obstacles in search- 
ing magic structures of group-IV nanostructures accord- 
ing to our model analysis. Firstly, we can estimate the 
Hamiltonian simply by E tot = m/ix +^Mh — (2mEx-x~ 
0.5nEx-x + nEx h), since the total energies (E tot ) di- 
rectly obtained from the first-principles calculations do 
not involve the environment-related chemical potentials 
and in fact correspond to the Hamiltonian in Eq. (1) 
with fix = Mh = 0- Secondly, we should only consider 
few possible candidates of X m H n with various n for a 
certain m, as isomeric structures with the same chemical 
formula will possess comparable total energies. Thirdly, 
we have searched the least n corresponding to a certain 
m iteratively and determined magic structures by the 
Simplex method. 

To verify the reliability of the above model, we in- 
vestigate the energetic stability of group-IV nanoctrys- 
tals (with the example of Dia-NCs and SiQDs) to 
search the magic structures, using the first -principles 
method implemented in Vienna Ab initio Simulation 
Package (VASP)0, [25|. We use Vanderbilt ultrasoft 
pseudopotentials[26] and the exchange correlation with 
the generalized gradient approximation given by Perdew 
and Wang [2 7]. We set the plane- wave cutoff energy to 
be 350 eV and the convergence of the force on each atom 
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FIG. 3: (Color online) Formation energies as a function of 
the number of Carbon atoms for diamond nanocrystals with 
/jlh = — 4eV . Inset shows the chemical phase diagram for 
C m H n with m < 12. 



to be less than 0.01 eV/A mesh of k space is used and 
the vacuum distance is set to be 9 A, which is enough to 
make the systems isolated. 

We define the cohesive energy (i£ C oh) per X atom in the 
nanocrystals X m H n as E coh = (m/i Xo +™Mh - E tot )/m, 
i.e., E coh = 2E X -x ~ (0.5£x-x - Ex-n)ot, with a is the 
H/X ratio (n/m). Figure 2 shows the cohesive energy of 
X m H n (X=C, Si) as a function of the H/X ratio. As pre- 
dicted, the cohesive energy decreases with the decrease of 
the H/X ratio following a linear relationship, approach- 
ing the value of bulk material ( — /iDiamond-buik = 

7.65eV 

and — /isi-buik=4.67eV[28|) when the H/X ratio reaches 
zero. In addition, we can make a deduction that the 
reaction heat Q = £ to t(C m H n ) - E to t(C m , Diamond) - 
£tot(H n ,H 2 ) = n(0.5£c-c - #c-h + 0.5£ H -h) will also 
have a linear dependence on the H atom number n, con- 
firmed by our calculations (shown in the inset of Fig. 
2). Thus, the interaction in X m H n is dominated by the 
localized X— X and X— H interactions and the Hamilto- 
nian estimated in the model by E tot = m/ix + n ^n — 
(2mEx-x — O.bnEx x + nEx yl) is at the accuracy level 
of the first-principles approach. 

In the following, we investigate the magic structures 
of Dia-NCs exemplified by C m H n . We calculated the 
total energies of nanocrystals with carbon atoms m < 
12 and obtained the formation energies as a function of 
hydrogen chemical potential by Ef = (E tot — m/ix — 
^Mh — nfiiij/m. As shown in the inset of Fig. 3, we find 
that there is a critical point //o(Mo = 0.5Eq-c ~ Eq-h — 
— 2.63eV) of chemical potential, at which the formation 
energies are the same for all these nanocrystals. Below 
the critical point C10H16 is the most stable, while CH4 is 
the most stable when /an is above the critical point. As 
predicted, CH4 and C10H16 are stable states when the 
number of carbon atoms m < 12. 

Figure 3 shows the formation energies for various 
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FIG. 4: (Color online) The possible dimer reconstructions and 
their effect on the number of hydrogen atoms.. 



C m H n with the chemical potential of hydrogen /in = 
— 4eV. The hollow circles are from our model prediction, 
which is in excellent agreement with the ones from the 
first-principles calculations (marked with solid triangles). 
It is noted that the isomeric structures have similar for- 
mation energies especially for the chemical potential far 
away from the critical point, though there are remarkable 
difference in their total energies. For a certain m, the for- 
mation energy decreases as the hydrogen atom number n 
decreases. The magic structures can be found at the lo- 
cal minimum of the formation energies, such as C10H16, 
C14H20, C18H24 and C22H28- All these structures are 
consistent with the previous model analysis, which have 
also been confirmed by experiments [!, We can also 
predict that, for the chemical potential /in > Mo, the for- 
mation energy will decreases as the hydrogen atom num- 
ber n increases for a certain m, thus the stable structures 
are alkane C m H 2m +2. 

For hydrogenated group-IV nanocrystals, magic struc- 
tures often correspond to the ones with the least or most 
hydrogen atoms, according to the linear dependence of 
total energies on H/X ratio. In our previous study [29], 
we have investigated magic structures of hydrogenated 
SiQDs, which are in agreement with the ones of Dia- 
NCs. It is reasonable as our model analysis showed that 
the magic structures of X m H n are the same for all group- 
IV elements. So far, we have not considered any possible 
reconstructions, which might induce strain and instabil- 
ity, especially for Dia-NCs with small size. For SiQDs, 
however, the formation of dimer s is common and often 
dominates the surface reconstruction [ll[ [l6|, [30|, which 
will further decrease the hydrogen number. According 
to our calculations, we find that the i? C oh of SiQDs with 
reconstructions also follow the linear dependence on the 
H/Si ratio. 

For simplicity, we consider a dimer reconstruction on 
(100) facet as is shown in Fig. 4. The hydrogen number 
of structure Si will decrease by 2 when the dimer forms. 
For structure S2, the hydrogen number will decrease by 
2 or 4, depending on the selection of Si atoms for recon- 
structions. To find out the max decrement of hydrogen 
number, we construct a matrix comprised of and 1 ac- 
cording to the arrangement of Si atoms with — SiH2i the 
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FIG. 5: (Color online) (Color online) Formation energies as a 
function of the number of Silicon atoms for SiQDs with and 
without the dimer reconstructions. 

element M(i,j) in the matrix is 1 when a dimer could 
be formed by the ith and jth Si atoms; otherwise, it 
is 0. Mi and M2 correspond to the structures Si and S2 
shown Fig. 4. The max decrement of hydrogen number is 
equal to the rank of the matrix. It can be comprehensible 
that the max decrement of hydrogen number is indepen- 
dent of the order of Si atoms, and the rank of matrix 
is an intrinsic parameter that conserves in the elemen- 
tary transformations. For SiQDs with various facets, we 
can construct a block matrix:M=(Mi,0,0;0,M2,0;0,0,...). 
The max decrement of hydrogen number should be cal- 
culated as:rank(M)=rank(Mi)+rank(M2)+... We search 
the least hydrogen number for SiQDs with a certain num- 



ber m through the similar procedure. 

Figure 5 shows the formation energies of Si quan- 
tum dots with and without dimmer reconstructions with 
/in = — 4eV. The dimer reconstruction decreases the for- 
mation energies and changes the magic structures. For 
example, SiioHi6 might not be the magic structure as its 
formation energy is higher than that of SisHi2, which can 
be obtained from SigHie with two dimer reconstructions. 
Besides, Sii 2 Hi 6 , Sii 5 Hi 8 , Si 2 oHi 8 , and Si 2 5H 2 o are new 
magic structures. The reconstructions also have signif- 
icant effect on the symmetry of SiQDs. Stable SiQDs 
without reconstructions tend to be octahedron enclosed 
by (111) facet, e.g. SissE^. However, the reconstructed 
Si2oHis has the symmetry of C$ v and Si25H2o is approx- 
imately spherical. 

In summary, we have proposed an effective method 
for investigating the stability of H-terminated group-IV 
nanostructures. We express the Hamiltonian of X m H n 
(X=C, Si, Ge and Sn) into the linear combination of the 
atom numbers and obtained the total energies with high 
reliability. We find that the isomeric structures are nearly 
energetically degenerated and the stable structures corre- 
spond to the ones with the least or most hydrogen atoms. 
The predicted stable Dia-NCs are in good agreement with 
experimental observations. With the overcome of the 
bottlenecks for stability determination, our model pro- 
vides an efficient technique of searching magic structures. 
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